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Abstract-This paper designs a hybrid genetic algorithm for job- 
shop scheduling combined with a Genetic Algorithm (GA), 
Simulated Annealing (SA), anda Tabu Search (TS) algorithm. It 
gives a method for initial population generation for avoiding 
unreasonable solutions, and introduces SA and TS into the 
genetic operation mechanism of GA for overcoming the lack of 
job-shop scheduling optimization. Simulated computations show 
the superiority of the hybrid genetic algorithm, and examples 
verify the algorithm's efficiency. 
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I. INTRODUCTION 
Because any job-shop scheduling algorithm has different 
degrees of advantages and disadvantages, researchers in 
recent years have combined various algori-thms to utilize 
advantages and to avoid disadvantages for achieving optimal 
scheduling. Typical algorithms include the choice of 
combining greedy and simulated annealing algorithms; an 
algorithm of this type searches by using a greedy algorithm, 
makes the obtained job as initial JIE, and then uses the 
simulated annealing algorithm to resolve a machine 
scheduling problem. Experimental results show that this 
method is better than using just simulated annealing or greedy 
algorithms [1], Other options combine simulated annealing 
with heuristic algorithms to resolve job-scheduling problem 
with delivery [2], and combine genetic with tabu algorithms 
to solve earliness/tardiness schedu-ling problem with waiting 
casts [3], Melt -Annealing algorithm, an improved simulated 
annealing method for flow shop problems proposed in 
reference [4] analyses asym-ptotic convergences to global 
optimal solution by six different random sampling patterns, 
and solves the sort of problem of Flow Shop with a minimal 
Makespan index, Unlimited Intermediate Storage (UIS), finite 
Intermediate Storage (LIS), and No Intermediate Storage 
(NlS).It is an improved GA method proposed in reference [5] 
by minimizing the maximum completion time, which 
introduces machine selection items in chromosome coding to 
improve equipment utilizations and balances, or the 
chromosome encoding method used in reference [6], which 
generates different codes according to various processes. 

Let numbers represent the process of one work piece. All 
these distinct numbers are coded for evolutionary computation, 
and a matrix decoding method is proposed using a multi- 
objective job-shop scheduling algorithm as in reference [7], 
which combines positive-reverse scheduling with a two- 
population genetic algorithm. 

To improve the local search ability of GA, also to improve 
further optimization quality and search efficiency, and 
compensate for the lack of a single optimization method, this 
paper gives a GA-SA-TS algorithm, which combines GA, and 



tabu search with simulated annealing algorithms. Examples 
are given to verify the efficiency of the algorithm. 

II. MATHEMATICAL MODE OF JOB-SHOP 

SCHEDULING 
Job-shop scheduling is a simplified model of many 
realistic production scheduling problems, and it is a typical 
NP-hard problem. Job shop is a production procedure using 
workshop resources for a certain object, while the job 
scheduling effectively sorts the job shop and minimizes (or 
maximizes) a certain objective function. Here, a job object 
can be a component, or a part of the component; it can be one 
part of a process or a certain process step. Because the process 
is a basic element of components and parts, the work of the 
same process step is usually an absolute division of process. 
That is, these will be finished on the same machine. Therefore, 
we established the job shop object as a part of the process. 
Certainly, job scheduling objects are established in the 
process as well. 

A. Basic Constraints of Job-shop Scheduling 

Job-shop scheduling not only has the extreme property, 
that is, its result is the minimum (or maximum) value of a 
solution, it also has the effectiveness property, that is, its 
result must meet the basic constraints of job-shop scheduling 
shown as follows: 

• Every process must be processed and produced on its 
specified machine. 

• Every process must wait for entering into a process 
procedure until its previous process has been finished. 

• A machine only can process one part at a time. 

• The procedure of processing parts on the machine is not 
allowed to be interrupted. 

B. Mathematical Model 

To describe the mathematical model of JSSP, we first give 
the following hypothesis and definition. 

Assume that P denotes a collection of n parts; then P can 
be described as P= {PI, P2, P3— , Pn}. 

Assume that M denotes a collection of m machines; then 
M can be described as M= {Ml, M2, •••, Mm}. 

Assume that the process number of Pi is Ki, and JPi 
denotes the process collection of Pi part; then JPi can be 
denoted as JPi= {Ji (1), Ji (2), ••• , Ji (k)}. 

Definition 1: Suppose the elements of process Ji belonging 
to Pion Mas:Mi(l),Mi(2)..., Mi(k), denoted by using i(l), 
i(2)..., i(k), where, k=max{kl, k2...,.kn}, 0^(i(j))^ m, 
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0— i — n, 0— j — k. Then, the matrix of a machine process 
order can be described as: 



Q= 



1(1) 1(2) 
2(1) 2(2) 

n(l) n(2) 



l(/c) 
2(fc) 

n(k) 



Definition 2: Suppose that the process order on machine 
Mi is JMi, where 0<i<m; then JMi can be denoted as JMi= 
{j(l),j(2)...t j (li)}. Here, li is the total number of process 
orders on Mi, and j(i) is the order number of the i' processed 
part. 

Definition 3: Matrix JM denotes all the arrangement of 
machine processes; we call JM to be the arrangement matrix 
of the machine process. Formally, it is: 



1(1) 1(2) 
2(1) 2(2) 

m(l) m(2) 



JM= 



l(wf 




JM,' 


2(w) 


or JM= 


JMi 


m(w) 




JMm 



Here, JMi denotes the process order on the i machine. 

Definition 4: We define the process time matrix T as 
follows: 



tl(l) 


tl(2) • 


■ tl(w) 


t2(l) 


t2(2) • 


• t2(w) 


tm(l) 


tm(2) ■ 


• tm(w) 



Each element in the matrix corresponds to an element in 
JM. Let (tj (w)) denote the finish time of j (w). 

Assume that P is the set of parts, M is the machine set, Q 
is the matrix of machine process order, and JM is the order 
matrix for the machine process. 

For a specified machine processes matrix JMO, if it meets 
the formula T(JMO) = minT(JM), that is, JMO minimizes the 
objective function F(JM), and compatibilities with Q, then we 
call JMO as JSSP optimal solution; call the procedure of 
solving JMO as job-shop scheduling. Its objective function is: 



F = min(T(J m )) = minmaX 



! max^ } 
i<y<v 



Here, w=max {kj}, k ; is the order number of workpieces J;, 
Tjj denotes the required process time of the j th process of the 
i workpiece. 

III. HYBRID GENETIC ALGORITHM 
As a kind of group optimization, a genetic algorithm 
optimizes the solution set through mutation or restructuring 
operations, and implements fast search operations ultimately. 
But, it has early-maturing and poor local search ability 
problems. Simulated annealing can be introduced into a cross- 
selection mechanism of a genetic algorithm to avoid the 



shortcomings of early-maturing, when high temperature (pre- 
GA) similar to the individual fitness of offspring produces 
similar crossover probability, and when temperature falls, the 
enhanced role of stretches, similar to the fitness of individual 
fitness differences, amplification allowing excellent 
individual advantage is more obvious [8]. Variation of 
pressure with TS to solve the problem GA, GA mutation 
probability is difficult to solve the contradiction selected, and 
because the tabu search algorithm to avoid the circuitous 
search capability, enhanced local search hybrid algorithm to 
improve the convergence speed. 

A. GA-SA-TS Hybrid Strategy 

GA-SA-TS hybrid strategy is a double-layer parallel 
search structure. At process level, the hybrid algorithm carries 
out GA, SA and TS search under various temperatures. It is a 
kind of double-layer serial structure. Here, the initial solution 
results of SA come from the cross evolution results of GA; the 
solution of SA obtained by Metropolis sampling process will 
become further evolution of the population of GA. The initial 
solutions of TS come from the evolution of variation of GA, 
the solution obtained by TS also become a future evolution of 
population of GA. At space laver, GA provides a parallel 
search structures, it transforms TS, SA into parallel TS, SA 
algorithms. Therefore, the hybrid algorithm introduced in this 
paper is optimizing populations from parallel structure. 

GA-SA-TS hybrid strategy of search behavior is 
controlled. It can be controlled by initial temperature of SA, 
anneal function, sampling frequency and the length of TS tabu 
table. Control of initial temperature, initial search behavior 
can control algorithm, controlling the temperature level can 
control the strength of the jumping ability of an algorithm. 
The jumping ability under high temperature is helpful for 
avoiding premature convergence. High temperature high jump 
is beneficial to avoid falling into local minima; chemotaxis 
optimization under low temperature is helpful to improve 
local search ability. Controlling the rate of temperature 
decreasing can control the decreasing amplitude of jumping 
ability and influence the smooth of search process. 
Controlling sampling frequency can control the search ability 
under various temperature levels, and influence the 
corresponding cross probability distribution. 

Controlling and setting the length of tabu table can 
improve the performance of the algorithm. In short, 
controllable performances enhance the low ability of the local 
search of GA, and overcome the premature shortcomings. 

B. Key Technologies and Main Steps of Hybrid Genetic 
Algorithm 

Using hybrid genetic algorithm to solve the key 
technologies of job-shop scheduling includes the designs of 
genetic operator and tabu search. 

1) Process Based Encoding Method: 

The encoding method need to be carefully designed when 
using genetic algorithm to resolve the job-shop scheduling. 
The adopted encoding method must ensure that the scheduling 
results are feasible scheduling solution. This paper selected a 
process based expression method after considering various 
factors involved in encoding procedure. This kind of 
expression method encodes scheduleing into a sequence of the 
process. It makes each gene to represent a process, each gene 
chain to represent a possibility of sorting a group tasks. 
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Assume that there are n processes in a group of tasks; it is 
not specified for each process from 1 to n in encoding process 
(that is, using different codes represents different process). 
Because of the limitation of processing technology, this pro- 
cess method results that chromosome produces infeasible 
scheduling solution, and increases difficulty in operation. 

Base on the reasons above, in this paper, we assign the 
same codes for the process belonging to the same part, and 
explain them according to the appearance order in given 
chromosome. For n workpieces and m machines problems, 
there are n=m genes in a chromosome; each workpiece 
appears m times in chromosome, each gene does not indicate 
a specific process of the workpiece, but indicate a process 
with upper and lower dependency. It is easy to observe that 
the arbitrary arrangement of a chromosome always produce a 
feasible scheduling. In the following, a practical example is 
given for explaining the coding ideas. 

Workpiece 1: (1, 1, 1, 2) (1, 2, 2, 3) (1, 3, 4, 4) (1, 4, 1, 6) 

Workpiece 2: (2, 1, 5, 7) (2, 2, 4, 8) (2, 3, 2, 3) (2, 4, 3, 9) 
(2, 5, 1, 4) 



Workpiece 3: (3, 1, 2, 
(3, 5, 4, 7) (3, 6, 1, 9) 



(3, 2, 4, 7) (3, 3, 1, 1) (3, 4, 5, 2) 



1) 



Workpiece 4: (4, 1, 3, 1) (4, 2, 2, 10) (4, 3, 1, 2) 
Workpiece 5: (5, 1, 2, 13) (5, 2, 3, 15) (5, 3, 4, 1) (5, 4, 3, 

Workpiece 6: (6, 1, 2, 3) (6, 2, 4, 3) (6, 3, 5, 6) 

Here, the symbols in brackets above denote the number of 
workpiece, the order of process of the workpiece, the machine 
tool number used by the process and the process time of the 
process respectively. 

A chromosome consists of a group of digits, using a digit 
denotes a process of the workpiece. However, the same digit 
which appears in different position in chromosome denotes a 
different process of the workpiece. So, the chromosome 
according to this example consists of a encoded string which 
includes four 1, five 2, six 3, three 4, four 5 and three 6, etc. 
such as 2213212514361235566543433, where the digits i 
appearing at j times denotes that the j' process of workpiece i. 
Therefore, all the process of ij can be denoted by using 
symbol i. As shown in the following table, we explain the 
chromosome in the example as a feasible scheduling. The 
chromosome in the instance can be explained as feasible 
scheduling in Table I: 

TABLE I 
CHROMOSOME EXPLANATION 



Gene 


2 


2 


1 


3 


2 


Process 


(2,1,5,7) 


(2,2,4,8) 


(1,1,1,2) 


(3,1,2,6) 


(1,3,4,4) 


Gene 


1 


2 


5 


1 


4 


Process 


(2,2,4,8) 


(1,4,1,6) 


(5,1,2,13) 


(1,3,4,4) 


(4,1,3,1) 


Gene 


3 


6 


1 


2 


3 


Process 


(3,2,4,7) 


(6,1,2,3) 


(1,4,1,6) 


(2,5,1,4) 


(3,3,1,1) 


Gene 


5 


5 


6 


6 


5 


Process 


(5,2,3,15) 


(5,3,4,1) 


(6,2,4,3) 


(6,3,5,6) 


(5,4,3,1) 


Gene 


4 


3 


4 


3 


3 


Process 


(4,2,2,10) 


(3,4,5,2) 


(4,3,1,2) 


(3,5,4,7) 


(3,6,1,9) 



2) Cross Operators as the Most Important Operators: 

For the chromosome individual based on ordering 
processes, it produces non-means encoding string if adopting 
traditional single-point crossover operators. That is, it 
produces infeasible scheduling solution. As the following, this 
paper adopts the crossover operators which combine single- 
point with order crossovers, and describe the procedure 
combining with practical examples given in the encoding 
procedure. 

Stepl. In the two parent chromosomes popchroml and 
popchrome2 (hereinafter, PI and P2), we randomly select an 
integer smaller than Number (total number of genes belonging 
to a chro mo-some), from the first gene to i' gene are consi- 
dered as crossover region, and the encoding string in this 
region is denoted with Tl and T2. Now, there are two parent 
chromosomes: 

PI: 5125323161333663144222545 

P2: 2553211612524445323363361 

Randomly select i=6, then the coding string of the 
crossover region is: 

Tl: 5 12 532 

T2: 2 5 5 3 2 1 

Step2. In PI, search all the gene values in T2 from left to 
right, set the gene value corresponding to the position (digits 
of italics add thickly) of PI to 0. Similarly, in P2, search all 
the gene values in Tl from left to right, set the gene value to 
as well. 

PI: 5100003061333663144222045 

P2: 2053201610024445003363361 

Step3. To concentrate the gene with zero value of pi and 
p2 into the crossover region by left shifting, while keep the 
other genes with a non-zero value the same relative order. 

PI: 0000005136133366314422245 

P2: 0000002532161244453363361 

Step4. To replace the content of PI crossover region with 
the content of T2, replace the content of P2 with the content 
of Tl. So far, two child individual childl and child2 can be 
obtained respectively. 

3) The Main Steps are Described as the Following when 
Applying Hybrid Genetic Algorithm to Solve Job-Shop 
Scheduling: 

Assume S is a solution set consisting of all possible 
solutions. For arbitrary i G S, it has a determined objective 
function value f (i). Schedule is used to save current optimal 
chromosome. TotalTime is an optimal individual fitness value. 

Here, the iterative target is solving minimum value of the 
function based on the feature of sche-duling problem. Its main 
steps are: 

Stepl. Set the size of initial population as pop-size, scale 
of crossover and variation with Pc and Pm respectively; set 
cooling schedule para-meters of SA under the concept of 
quasi-equilibrium: initial temperature TO, temperature 
attenuation coefficient a, set the number of can-didate 
solutions and tabu table length of TS algorithm dynamically, 
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set the adaptive optimal scheduling best and coincided 
Fit Value. 

Step2. Using for process based decimal encoding products 
popSize initial feasible solutions, let the initial search algebra 
gen = 0. 

Step3. Application of heuristic rules on the chromosome 
decoding, arranging the process order of each part to be 
processed on each machine, computing and evaluating the 
fitness function value of the current population. 

Step4. Determine whether a termination condition is 
satisfied. If the termination condition is satisfied then the 
algorithm ends else continue. 

Step5. Let Tk + i=aTk, gen + =l, according to crossover 
probability Pc, using simulated annealing crossover operators 
products new individuals, to decode and compute fitness 
value for it. 

Step6. On the solution to the mutation probability Pm for 
screening, for each individual tabu mutation screening: 

6-1. fromthe current population selected solu-tion into the 
tabu search stack. 

6-2. pop a solution x[i] from tabu search stack, set tabu 
table as null, optimal state as null. 

6-3. dynamic generation neighbourhood solution and the 
number of candidate solutions from the current solution. 

Check whether candidate solutions meet contempt criteria 
or not. If meet, then replace the optimal state y which meets 
the contempt criteria with x[i] , forms current new solutions, 
that \s,x[i]=y, put y corresponding tabu target into tabu table, 
and release the targets which meet tabu length len from tabu 
table, replace optimal state with y, then return to step6, 
otherwise return to step5. 

6-4. check candidate solution corresponding to tabu 
features; select the non-tabu target in candidate solution set 
corresponding optimal state as new current solutions. 
Meanwhile, put the correspondding targets into tabu table, 
and release those targets meeting tabu length fromtabu table. 



6-5. if the required number of tabu search has been met, 
then ends the searching for x[i], jump to step7, otherwise 
return to step3. 

6-6. if all the solutions in the tabu search set have been 
carried out in tabu search, then stop carrying out the tabu 
search algorithm, otherwise return to step2. 

Step7. Save the current solution with best performance 
into Best, fitness value into Fit Value. 

Step8. Return to Step4. 

Explanationl. The number (candi) of candidate solutions 
in this paper is related to the scale of the problem, it is 
generated dynamically. The following is the rules: 



candi = 



n 

In 

lOn 



0<n<25 
25<n<50 
50 < n < 100 



4 


< n < 20 


5 


20 < n < 60 


6 


60 < n < 100 


7 


n =100 



Explanation2. Dynamic generation tabu table length len as 
the following rules: 



len =i 



C. Experiments 

For small scale production like data MT10, the sorted 
GANTT figure by device sequences can be obtained after 
system running, as shown in Fig 1. 

To investigate the performance of hybrid strategy of GA- 
SA-TS, and compare GA, SA, TS, we simulated 11 typical 
examples for 30 times through example simulation module. 
Table II summaries the obtained optimal solutions and 30 
times simulated average values. 
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Fig. 1 Sorted GANTT figure by device sequences 
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TABLE II 
Performance comparison ofGA-SA-TS, GA, SAandTS 



Problem 


C** 


GA-SA-TS 


GA 


SA 


TS 








C* 


% 


C* 


% 


C* 


% 


C* 


% 


MT06 


55 


55 





55 





55 





55 


3.273 


LA01 


666 


666 





666 





666 





666 


3.048 


LA06 


926 


926 





926 





926 





926 





LA11 


1222 


1222 





1222 





1222 





1222 





MT10 


930 


930 


2.34 


997 


11.87 


939 


6.65 


962 


10.89 


MT20 


1165 


1165 


1.12 


1247 


12.40 


1227 


8.33 


1238 


11.41 


LA16 


945 


945 


0.5 


979 


6.032 


979 


4.17 


979 


5.032 


LA21 


1046 


1046 


1.147 


1156 


14.56 


1083 


7.49 


1157 


15.53 


LA26 


1218 


1220 


1.181 


1328 


12.78 


1253 


5.10 


1290 


10.78 


LA31 


1784 


1790 


0.785 


1836 


4.826 


1784 


0.39 


1785 


1.234 


LA36 


1268 


1268 





1384 


11.91 


1321 


5.80 


1302 


4.236 



The 11 typical examples include: MT06(6><6), 
MTlO(lOxlO), MT20(20x5) and LAO 1(10x5), LA06(15x5), 
LAll(20x5), LA 16(10x10), LA21(15xlO), LA26(20xlO), 
LA31(30xlO), LA36(15xl5), where MT and LA represent 
question constructors Mute, Thompson and Lawrence, and the 
digits in brackets represent the number of parts and the 
number of machines respectively[9]. 

On the optimization of GA-SA-TS indexand GA, SA and 
TS for comparison, the statistical results are shown in Table II. 
In the Table, the meaning of each symbol, C * * is for the 
optimal solution of the problem. The C * algorithm simulation 
of the optimal value is of 30% average optimal value relative 
to C * * deviation. 

IV. CONCLUSIONS 

In this paper, the researched hybrid genetic algorithm 
selected practical problems as target; it has a good effect on 
solving multi-workpiece, multi-process and multi-machine 
mixing flow processes; it also has a better application value 
for improving the efficiency of job-shop scheduling. 
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